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ABSTRACT 

y 

Optical  fibers  are  produced  by  heating  a  silica  rod  and  pulling  on  the 
rod  to  produce  a  thin  fiber.  This  process  involves  several  interacting 
effects,  especially  the  temperature  dependent  viscosity  of  the  rod,  and  the 
shape  of  the  rod  in  the  transition  zone.  In  this  paper  the  equations 
describing  this  process  are  studied  by  means  of  a  perturbation  analysis  using 
as  the  perturbation  parameter  the  ratio  of  the  rod  to  the  length  of  the 
furnace.  The  silica  rod  is  usually  doped  so  that  it  will  have  special  optical 
properties  and  it  is  important  to  know  the  concentration  of  the  dopant  in  the 
final  fiber.  Simple  equations  for  the  significant  terms  in  the  perturbation 
series  are  obtained  for  all  quantities  of  interest.  A  simple  formula  is 
obtained  which  expresses  the  radius  of  the  fiber  in  terms  of  the  initial 
radius  of  the  rod,  the  drawing  force  and  other  important  parameters. 
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A  STUDY  OP  THE  DRAWING  PROCESS  FOR  OPTICAL  FIBERS 


John  C.  Strikwerda 


1.  Introduction 

Optical  fibers  are  usually  produced  f roe  silica  rods  by  heating  and  pulling  the  rod 
in  the  axial  direction*  This  essentially  sinple  physical  process  is  difficult  to  analyse 
eathenatically  due  to  the  interaction  of  several  effects.  These  effects  include  the 
teeperature  of  the  furnace,  the  tasparature  dependent  viscosity  of  the  rod,  the  drawing  or 
pulling  force  on  the  rod,  and  the  shape  of  the  rod  during  the  drawing  process.  In 
addition  the  silica  rod  is  usually  doped  so  that  it  will  have  special  optical  properties 
and  it  is  important  to  know  the  radial  dependence  of  the  dopant  in  the  final  fiber.  In 
this  report  this  process  is  studied  so  as  to  produce  convenient  equations  fro*  which  to 
calculate  many  of  these  quantities  and  to  demonstrate  the  essential  features  of  this 
process.  In  a  subsequent  report  these  equations  will  be  solved  numerically  to  produce 
quantitative  results. 


2.  Description  of  the  Problee 

Consider  a  circular  rod  placed  in  a  furnace,  which  we  assume  as  cylindrical  in  shape, 
with  teeperature  distribution  in  the  furnace  given  by  ( Z ) .  The  rod  and  furnace  are 

taken  as  coaxial  with  the  Z-axis ,  and  for  simplicity  of  analysis  both  the  rod  and  furnace 
are  assumed  to  be  infinite  in  extent.  (In  general,  the  furnace  will  be  cold  for  all  but  a 
finite  interval.)  The  rod  is  forced  to  move  in  the  axial  direction  under  the  influence  of 
a  given  drawing  force  Fd.  The  initial  feed  rate  of  the  rod  is  U0  and  the  radius  is 
Rq.  Because  of  the  temperatures  of  the  furnace  the  rod  will  behave  as  a  viscous  liquid 
and  be  drawn  down  to  a  radius  Rj  and  have  velocity  U.,.  An  objective  of  this  study  is 
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to  show  how  the  temperature  distribution  Tf  (Z) ,  the  drawing  force  Fd,  and  feed  rate 
UQ  affect  the  radius  reduction  from  RQ  to  Rj.  The  rod  and  furnace  are  assumed  to  have 
no  angular  variation,  which  means  all  effects  will  depend  only  on  the  axial  and  radial 
coordinates.  Paek  and  Runk  ( 1978)  presented  a  method  for  determining  the  neck-down  shape 
of  the  rod,  however  there  are  several  difficulties  with  the  equations  which  they 
employed.  The  equations  should  allow  for  the  degenerate  case  in  which  the  rod  and  furnace 
are  at  the  same  constant  temperature.  The  equations  of  Paek  and  Runk  do  not  have  such  a 
solution  unless  the  temperature  is  1100*C,  which  appears  to  have  been  chosen  quite 
arbitrarily.  The  equations  used  here  do  not  have  such  arbitrariness.  Because  of  the 
greater  generality  of  the  equations  employed  here  they  should  be  useful  in  studying  the 
cooling  of  the  fiber  after  it  exits  from  the  furnace.  The  cooling  process  is  of  interest 
because  of  its  role  in  determining  the  strength  of  the  fiber. 

The  method  used  in  this  analysis  is  the  method  of  perturbation  expansions  using  the 
ratio  of  the  rod  radius  to  a  length  representative  of  the  heating  zone.  This  Blenderness 
ratio  is  of  the  order  of  0.05  for  the  draw-down  process. 

The  use  of  the  slenderness  ratio  has  been  employed  to  study  slender  jets  of  inviscid 
fluids  with  great  success,  see  the  work  of  Geer  ( 1977a, 1977b)  and  has  led  to  an  effective 
numerical  method,  see  Strikwerda  and  Geer  (1980).  (See  also  Geer  and  Strikwerda  (1980)  and 
(1983).)  An  advantage  of  the  perturbation  method  is  that  it  displays  the  relative 
importance  of  the  various  effects  which  interact  in  the  process.  For  the  draw-down 
process  the  perturbation  method  gives  relatively  simple  formulas  for  the  determination  of 
the  temperature,  shape,  and  dopant  concentration.  An  especially  striking  formula  is  that 
for  the  rod  shape  as  a  function  of  the  viscosity  of  the  rod. 

3.  The  Governing  Equations 

The  governing  equations  will  be  given  first  in  dimensional  form  and  then  rewritten  in 
nondimensional  form.  The  convention  which  is  adopted  is  that  for  quantities  denoted  by 
Latin  letters  the  dimensional  quantity  will  be  denoted  by  the  upper  case  and  the 
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nondimensional  quantity  by  the  lower  case.  The  exception  is  temperature  for  which  0 
will  denote  the  nondimensional  temperature.  For  quantities  denoted  by  Greek  letters  an 
overbar  will  distinquish  the  dimensional  value.  The  exception  here  is  that  e  denotes 
the  emissivity  but  e  is  the  slenderness  ratio  of  the  perturbation  expansion. 

The  equations  describing  the  drawing  process  are  essentially  those  given  by  Paek  and 
Runk  (1978).  The  differences  between  the  equations  sets  are  that  in  this  paper  the 
pressure  is  included  in  the  stress,  the  heat  flux  is  calculated  differently,  and  here  the 
dopant  concentration  equation  is  also  included. 

The  equations  governing  the  movement  of  the  rod  are  the  balance  of  forces 


3  o  3(Rcr  ) 

zz  _  rr 

3Z  R  3R 


3o  .  3<R a  )  a,. 

rz  _1  _  rr  td 

3z  R  3R  R 


where  the  a ^  are  the  usual  coaiponents  of  the  stress  tensor.  At  the  surface  of  the  rod 


we  have  the  conditions  of  no  stress 


a  n  +  a  n 
it  t  rs  r 


°rz  "*  +  ®rr  nr  ’  0 


where  (nE,nr)  are  the  coordinates  of  the  unit  normal  at  the  surface.  Surface  tension  is 
ignored  as  it  is  negligible  for  silica  rods  (Geyling  and  Homsy  1980)). 


The  rod  is  an  incompressible  medium  so  the  velocity  components  satisfy 

3D  1  >R» 


—  +  -  m  0 

3z  R  3R 


The  liquid  phase  is  assumed  to  be  a  Newtonian  fluid  (see  Manfre  (1969))  so  the  stress 
components  in  the  liquid  are  given  by 


—  -  3u 

-  P  +  2W  3Z 


-  -30  ^  3v, 
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where  P  is  the  pressure  and  u  is  the  viscosity.  Gravitational  effects  can  be  included 
in  the  pressure,  or  ignored  as  we  do  here  (Geyling  and  Homsy  (1980)).  The  equation  for 
the  temperature  is 


in  the  rod  and 

on  the  surface,  where 
Q(Z)  is  given  by 


4n2  o 

P  C  (0  ♦  v  )  -  V2(K  T  -t  — —  T4) 

v  3Z  3r  c  _  — 

3  Y 

,  4n2  a 

s-  (K  T  ♦  -----  T  )  -  -  Q(Z) 

3 7 

Q(Z)  is  the  heat  flux  due  to  the  furnace  ( Paek  and  Runk  (1978)). 


_ •  >  .  .  cos*  cos* 

Q(Z)  -  e  a  I  ]  °  (T(Z)  -  T  <t)  )  - r - -  Ad*d£ 

-  -*0(z.C)  *  R'2 


where  R12  is  the  distance  between  the  point  P1  on  the  rod  surface  at  Z  and  the 
point  P2  on  the  furnace  surface  at  £,  (see  Sparrow  and  Cess  (1967)).  The  angle  *1 
is  that  between  the  normal  at  P.,  and  the  line  connecting  P,  to  P2,  and  +2  is  the 
angle  between  that  line  and  the  normal  at  P2.  The  radius  of  the  furnace  is  A  and 
*q(Z,C)  is  the  angle  intercepted  by  the  tangent  plane  at  P1  on  the  section  of  the 
furnace  at  ? .  The  other  quantitites  in  the  above  equations  are  the  density  p,  the 

specific  heat  Cy,  the  thermal  conductivity  Kc,  the  index  of  refraction  nQ,  the 
absorption  coefficient  T,  the  esdsaivity  e,  and  the  Stefan-Boltsnann  constant  a. 

The  equation  governing  the  transport  of  dopant  is 

DH  +  v  fi  ■  k  (K<T>  If*  (RK(T)  fi> 

where  C(R,Z)  is  the  concentration  of  dopant  and  K(T)  is  the  diffusivity.  On  the 
surface  the  condition  is 

The  functions  that  are  presumed  known  are  the  viscosity  u (T)  and  diffusivity  K(T), 
as  functions  of  the  temperature ,  and  the  temperature  of  the  furnace  wall  Tf(z)  as  a 
function  of  the  axial  distance  Z. 
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The  variables  ar«  nondimensionalised  by  expressing  then  as  multiples  of  various 


constant  reference  quantities.  The  reference  velocity  is  the  feed  rate  UQ,  the  reference 

force  is  the  drawing  force  Fd.  the  reference  temperature  is  the  initial  temperature  of 

the  rod  Tq  in  an  absolute  scale ,  e.g.  *K.  The  initial  furnace  temperature  is  assumed  to 

be  the  same  as  that  of  the  rod.  There  are  two  reference  lengths,  one  being  the  initial 

tod  diameter  RQ,  the  other  being  a  representative  length  of  the  furnace  ZQ.  The  ratio 

I^/Z0  is  assumed  to  be  sswll. 

The  non dimensional  variables  are  defined  by 

c  -  R/Z  r  ■  R/R  s  -  Z/Z 

o  o  o  o 

u  ■  o/o  v  ”  V/(0  t) 

o  o 

0  -  T/T  k  (ft)  -  K(T)Z  /OR2 

o  o  o  o  o 

The  reference  stress  is 


and  so 


o  -  f./r 
o  d  o 


®ij  “  Oij/Oo 

p  -  U  0/(0  Z  ) 
o  o 

p  ■  v/a  . 


The  concentration  is  nonalised  so  that  initially 

A 


j  C(R,Z)RdR  -  V2A  C 


and  c  «  C/C0. 


The  non  dimensional  equations  which  hold  in  the  rod  are 
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In  the  liquid  phase  of  th«  rod  tha  straaaas  are  given  by 

(6)  “  -  P  ♦  2M  £ 

-1  .3u  2  3v, 

(7>  °r«  “  e  »,<37  +  e  3l> 

(8)  °rr  “  '  p  +  2M  fe 

(9)  ■  -  p  ♦  2|i  v/r. 

The  normal  to  the  rod  surface  is  given  by 


/lecVc*)2  '  r 


where  the  rod  surface  is  given  by  r  •  s(s). 
The  boundary  conditions  are  therefore 


®rs  *  *  •*(*)  °s. 

«rr  “  e  ••<«> 

V  -  s’(s>  U. 


(fj  -  c2s'  |j>  (k1  ©  ♦  k2  o4)  -  -  q(sH1+€2s*(s)2)^2 


-  ♦  (s,5) 


cos*. cos*. 


q(s)  -  c  k,  a  j  J  °  <©(s)4  -  ©.(C)4)  - \ - 1  «♦<*■ 

-  -*.<8,0  v  r2 


<«>  ff-'2*’  If 

Because  of  the  nondiaensionalisation  Which  is  eeployed,  the  values  of  the  several 
variables  in  tha  initial  solid  phase  of  the  rod  (i.e.  at  s  -  -*•>  ares 
u(r,-»)  -  1f  v(r,-*»)  *0,  0  (r,^*>  ■  1 

|i(1)  -  »,  ke<1)  -  0,  axx{r~m)  ”  1 


and  tha  other  stress  coeponents  vanish 


and  we  assume  that  q(z)  vanishes  as  z  ♦ 


Also  c(r,-**)  is  some  specified 


function .  The  initial  radius  of  the  rod  is  1. 

If  0*  is  the  nondimensional  softening  temperature  of  the  rod,  then 


if  0  <  0* 


k  (0)  -  0 
o 


U(6)  <  • 


k  (9)  >  0 
o 


if  0  >0* 


A  possible  expression  for  u<0)  for  0  >  0*  is  u(9>  “  u*eb/^9”9  * '  Bee  Napolitano  and 


Hawkins  (1964).  The  nondimensional  constants  are 


k.  -  (KZ)/(pC  0  R*) 
1  CO  ▼  o  o 


k  -  (4?  3  n2  OZ  )/<3YPC  DR2) 
2  o  o  o  v  o  o 


k  -  (60  T  3)/(pC  D  Z  ) 
3  O  V  o  o 

a  -  A/R 


4.  Perturbation  Analysis  for  the  Shape 

Sach  of  the  variables  in  the  process  are  expanded  in  perturbation  series. 


For  example , 


u(r,z)  -  u  (r,z)  ♦  e  u,(r,z)  +  c  u,(r,z>  + 
O  1  2 

s(z)  -  s  (*)  +  e  s  (a)  ♦  e2s  (z)  ♦  •••  , 

0  12 
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Bimuh  of  the  subscripts  on  the  stresses  the  perturbation  index  for  the  stresses  will  be 


a  superscript,  and  a ^  will  be  the  only  variable  with  a  index  of  -1,  i.e. 
-1-1  o  o 

a  -  e  c  +eo  +•••• 
rs  rs  rz 


From  equation  ( 1 )  it  follows  that 


3ro 


r3r 


3ra 

rz 

r3r 


-  0, 


thus 

that 


ro 


-1 

3 

rz 

-1 


and 

and 


T<>rt  are  independent  of  r . 


The  boundary  condition  (10)  further  shows 


vanish  identically. 

-1 

a  • 
rs 


Since 

3u 
_ o 

uo  3r 


-  0 


it  follows  that  uQ  is  a  function  only  of  z.  From  the  continuity  equation  (3)  and  (12) 
there  follows 


du  (z)  s2 
o  o 

dz  2 
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r  dr  • 


s  3rv 
i  o  _ o  . 

L  —*r 


■  -  s  (z)v  -  -  s  ( z)  s'(«)  u  ( Z)  • 
o  o  o  o  o 

Therefore,  since  uQ(z)  -  1  in  the  initial  solid  phase, 

(17)  Uq(z)  -  1/s0(z)2  . 

Frost  (3)  and  (12)  it  follows  that 

(18)  v0(r,z)  ■  r  s£(z)/s0(z)3. 

Bquation  (1)  integrated  over  a  section,  together  with  (10)  yields 


s  (z)  3o 

J  °  — rdr  - 
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which  can  be  rewritten  as 
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.  s  (*) 

~  (J  °  a  °  rdr)  ■  0. 

d*  J  o  “ 

Baaed  on  the  initial  solid  phase,  it  follows  that 


s  <s) 

I  °  -1/2 

0 


Fro*  aquation  (2),  with  a ^  -  0, 


*  o  © 

I  Iflu  .-ii.  o 

rdr  r 


end  substituting  (8)  and  (9)  gives 


s—  ( -p  +  2  —  u  )  ■  t—  o  "0. 

dr  *o  .  .3  o  dr  rr 

s  (s) 


Since  (11)  shows  that  o  vanishes  on  the  surface  of  the  rod,  the  conclusion  is  that 

rr 

o  . 

a  m  o 

rr 


p  (r,s)  -  2 - ,  u  (r,s). 

O  .  v3  O 

■o(,) 


By  the  definition  of  o  in  the  liquid  phase,  with  (17)  and  (20) 


o  o 

o  rn  -  p  ♦  y  T— ^ 
ZS  rO  O  dz 


-  6  y  (r#z)  ■•(«)/•(«)" 
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At  locations  at  which  the  rod  is  entirsly  in  the  liquid  phase,  this  formula  with  (19) 


3  •o(*> 

-  6  s£(s)/so(«)  I  U0<r,s)  rdr  “Vj* 
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Ir*  l.v-*  *.• 


j'rs'ri 


In  the  intergral  In  r  make  tha  aubatitution 

(21)  r  -  a  (a)  n, 

o 

which  will  provi  wry  useful  liter  as  wsll,  to  obtain 

(22)  77^  -  -  —  (  v0(n,s)ndfi)  . 

At  this  point  it  ia  convenient  to  nake  an  approximation  by  raplacing 

p  (r,i)  by  u  (a),  which  ia  independent  of  r.  Tha  valua  of  |i  (a)  can  ba  obtained  in 
o  o  o 

any  of  aevaral  way a.  With  thia  approximation  (22)  can  ba  aimply  integrated  to  give  tha 
raaurkable  formula 

(23)  ”  **p  *“  6  ^  ME)'1  dc). 

It  ia  important  to  point  ant  that  tha  validity  of  tha  approximation  of  p  by  p 

o 

can  ba  checked  by  the  reenlta  of  oamputationa .  If  the  computations  ahow  that  for  a 

appreciable  interval  of  a  valuea  there  are  aubatantial  aolid  and  liquid  portiona  for 

each  value  of  a  then  the  above  approximation  ia  not  valid.  In  that  case,  if  the  rod  ia 

aolid  for  r  leaa  than  a  (a)  a  (a),  equations  (22)  is  replaced  by 

o  o 

t 

•  1  1  1 
(2«)  - 2 - -  --■*-</  v  (n»s)n«n) 

a  (1-(«  (a)a  (a))  )  a  (a) 

o  o  o  o 

which  follows  directly  from  (19)  and  o  °  equal  to  1.0  in  the  solid  phase. 


5.  The  effect  of  Dram-down  on  tha  Dopant  Concentration 
The  equation  governing  cQ(r,a)  ia 


It  is  advantageous  to  change  coordinates  from  (r  ,*)  to  (n<z)  where  n  is  defined  by 
(21)*  The  resulting  equation  is 


)c  -  .  9c 

ir"  "  n  9n  <n  W  liT 


z  >  0  <  n  < 


which  is  a  simple  parabolic  differential  equation  for  cQ.  The  initial  conditions  are 
that  co<n.»>  is  specified  in  the  solid  portion  of  the  rod,  (i.e.  where  n  “  r).  The 
boundary  condition  (15)  becomes 

9c 

(26)  -  0  at  n  -  1. 

Thus,  the  change  of  c0  is  not  explicitly  dependent  on  the  rod  shape,  depending  only  on 

the  temperature  9  ,  which  does  depend  on  the  shape.  Equations  (25)  and  (26)  are  very 
o 

amenable  to  numerical  computation  since  the  n  variable  la  on  a  fixed  domain. 


6.  Determination  of  the  Temperature . 

The  equation  for  the  nondimensional  temperature  0^(r,s)  is 

do  do 

u  ?— ™  ♦  v  ■  —  t—  (r  f—  (fc.0  +h 0*)) 
o  9s  o  9r  r  9r  9r  1  o  2  o 


Using  the  transformation  (21)  as  in  the  previous  section  this  becomes 

90  .  , 

<27>  SI*  -  ife  SI  «k10o+k2eo» 


The  boundary  condition  ( 13)*  becomes 


The  determination  of  q^z)  requires  some  analysis  to  obtain  the  limit  as  e 
saro. 


For  convenience,  set 


F(z,C,e)  -  04<sU),z)  -  e4(C> 


then  from  ( 14) 
(29) 


-  * 


cos*  cos*. 


q  (s)  -  k  a  lim  e  J  J  °  P(z,C>  - “ 

e+0  -•  -4  *  r 


2 

12 


d*dC. 


The  coordinates  of  the  point  on  the  rod  surface  at  s  can  be  taken  as  (: 

in  cartesian  coordinates  and  those  of  the  point  P2  on  the  furnace  surface  at 
(C,ea  cos  +,e  a  sin  4><  Thus 

r^2  *  (z-C)2  ♦  e2(s“a  cos  4) 2  +  c2(a  sin  4) 2 
«  (z-C)2  ♦  e2(s2  ♦  a2  -  2  as  cos  4) 

The  normal  to  the  rod  surface  at  P1  is 

(-es'(s),  1,  0)/(U(es*(*))2)1/2 
and  that  to  the  furnace  surface  at  F1  is 

(0,  -  cos  4,  -  sin  4) 

The  unit  vector  in  the  direction  P^j  is 

(C-*.  c(a  cos  4  -  s>,  e  a  sin  4>/*12 


Therefore 


cos 


*  _  £  <(s-C)s'(s)  ♦  a  cos4  -  s(z)) 
1  r12  / _ 2 


1+(ES' (Z) ) 


cos  4«  ■  — 
2  r 


(a  -  s(z)  cos  4>> 


12 


The  quantity  for  which  the  limit  is  taken  in  (29)  is  first  transformed  by 


tends  to 


,es(z),0) 
C  are 


replacing 


C  by  z-C,  and  for  convenience  set 


r  -  (a2*a<a)2  -  2as(s)  coe*)^2 


which  glvaa 


ks  a  c  •  *o  r(i,i-0(c«'+a  coa*-a)(a-a  coa*) 

^77a  L  L  ,,».«»  ?,» 

o 


»  <H<«a*>  )' 


ira  a  and  a'  ara  avaloatad  at  i< 
Naxt  aaka  tha  aobatitution 

tha  abora  quantity  la  than 


C  -  C  *  r. 


k,  a  •  f 


3  j  jo  F(a,a-«gCncga'g4a  coa*-a)(a-a  coa jj  ^ 

*(H(i»t7j  —  -d  r3  (C2v1>2 


Mow  tha  lladt  can  ba  takan  aaaily  giving 


"a  Mg 

c  -  “  r  <a,a)  J  — ,  ,  J  °  “r — — 

—  (?2-M)2  -♦  la,a>  (a2+a2- 


*  (a, a)  (a  coa* -a  ) (a-a  ooa* ) 

5 - ad*. 


*) 


3/2 


tha  firat  intagral  la  avaloatad  aaaily  aa 


_<s _ 

(C2+1)2 


-  «/2 


oaing  for  axaapla  C  *  tan  a.  tha  aacond  intagral  la  alao  aaaily  avaloatad  aa  follova. 

Firat  nota  that  *  (a, a)  la  givan  by 
o 

•  *"•  ♦«  ■  V8> 

alao  aatting  0  ■  a  (a)/a  tha  Intagral  baooaMB 
o 
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2  tco*0-8)(1-8oos»)  M  _  - 2  sln_j  ^  J4o 

0  (1+02-20cosO)3/2  (H-B2-2Bcos*/2  '° 

2  sin  ♦o  2  /T? 

"  *  2  • 

( 1+0  -2BcosO  ) 2  ( 1+0  -2B  )  2 

o 

therefor* 

(30)  qe  -  kj  ro  (*,*> 

■  k,  (04(s  (s),s)  -  Q*(z)) 

J  o  o  t 

Using  the  transformation  (21)  the  boundary  condition  (27)  is 

<31>  fe«kiWJ  *  -■ou>k3‘GJ<1**)  • 

aquation  (27)  with  boundary  condition  (31)  is  again  a  parabolic  initial  boundary 

value  problem  for  the  temperature  0  in  terms  of  the  coordinates  (n.x)  with  a  being 

o 

the  time-like  variable,  the  temperature  depends  on  the  fluid  flew  through  the  shape 
function  s0(a)  which  appears  in  (31)  and  is  determined  by  either  (22)  or  (24)  which 
depend  on  temperature  through  the  viscosity,  the  validity  of  the  use  of  an  average 
viscosity  |i(s)  Should  be  apparent  from  the  temperature  distribution  of  the  rod.  If 
there  are  significant  variations  in  temperature  in  the  rod  when  the  surface  begins  to 
bsuosm  liquid  then  the  use  of  equation  (24)  would  be  more  appropriate  than  (22). 

In  interesting  feature  of  the  formula  (30)  is  that  to  within  O(e)  the  flux  is 
independent  of  the  furnace  radius  a.  It  is  consistent  with  one's  intuition,  however, 
that  the  effective  furnace  temperature  is  approximately  that  of  the  furnace  closest  to  the 
rod.  However,  on*  might  expect  that  (30)  would  be  a  less  accurate  approximation  than  it 
is. 


this  section  serves  as  a  summary  of  the  previous  analysis.  The  primary  problem  is  to 

determine  the  temperature  distribution  9  (r.s).  The  governing  equation  is  (27) 

99  ° 

jf "  i  h fe  (kiVk2ei)) 
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for 


<*<•»,  0<r\<1 


with  the  boundary  condition  (31)  on  n  “  1# 


h  ■  -•o(*)k3(ei(1'*>  - 


The  surface  shape  function  s0(t)  is  given  by  (24). 


- 2 - 2  -  -  {  (  j  u(0  (n.s))  ndn)"1 

s  (x)(1-(a  (z)s  (z))  )  a  (s)  c 

o  o  o  o 


where  a  (z)  is  that  value  of  n  for  which 
o 


9  <  9*  for  n  <  a  (z) 
o  o 


0  >0*  for  n  >  a  (z) < 

o  o 


If  one  uses  an  approximate  viscosity  u(z)  then  (23)  holds 

£  m 

s (z)  ”  exp( -  ^  J  v(C)  ’dC). 


In  the  numerical  computation  of  these  problems,  one  would  take  the  furnace  temperature  and 


initial  rod  temperature  as  equal,  for,  say,  s  <  0,  i.e. 


0  (n,z)  -  0,(s)  ■  1  for  *  <  0. 

O  I 


Simultaneously,  the  dopant  concentration  can  be  computed  by  equation  (25)  with  boundary 


condition  (26) 


~  -  Jfc  VV"'*»  0<,“  1 


"0  at  n  ■  1. 

It  is  also  of  interest  to  obtain  expressions  for  the  drawing  radius  ratio  in  terms  of 


the  dimensional  parameters.  If  a  global  average  viscosity  u*  is  defined  by 


(U*)"1  -  ~  j  U(T(8))_1 


then  (23)  yields 


(32)  exp  <-Vo/<»oRo‘,#))* 


P 


k 


mmmmmmmwwm 


this  forauls  clssrly  show  ths  relationship  that  holds  between  the  drawing  force  Fd,  the 
feed  velocity  UQ,  the  initial  radius  RQ,  the  furnace  length  ZQ,  the  average  viscosity 
1 1*,  and  the  radius  ratio  R^/RQ.  Although  u*  depends  on  F^,  UQ ,  and  Rq,  presumably 
it  is  moat  strongly  dependent  on  the  furnace  temperature  distribution.  Assuming  this, 
formula  (32)  can  be  used  to  estimate  the  change  in  R1  which  would  result  from  small 
changes  in  F^,  Rq,  vo,  and  the  furnace  temperature. 

the  advantage  of  the  perturbation  analysis  is  that  it  gives  such  simple  expressions 
as  (32)  between  the  large  number  of  parameters  which  govern  this  problem.  The  relation 
(32)  is  only  true  to  within  terms  of  order  R0/ZQ,  but  the  experimentally  obtained  values 
for  the  coefficients  and  functions  which  have  been  assumed  as  known  may  well  have  errors 
of  a  similar  magnitude. 
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